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Critical flux pinning and enhanced upper-critical-field in magnesium diboride films 
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We have conducted pulsed transport measurements on c-axis oriented magnesium diboride films 
over the entire relevant ranges of magnetic field < iif < Hc2 (where Hc2 is the upper critical field) 



and current density 



< 1 < 



jd (where jd is the depairing current density) . The intrinsic disorder of 



the films combined with the large coherence length and three-dimensionality, compared to cuprate 
superconductors, results in a six-fold enhancement of Hc2 and raises the depinning current density 
jc to within an order of magnitude of jd- The current-voltage response is highly non-linear at all 
fields, resulting from a combination of depinning and pair-breaking, and has no trace of an Ohmic 
free-flux-flow regime. 



PACS numbers: 74.25. Sv, 74.25. Fy, 74.25.Bt 



I. INTRODUCTION 



gnesium diboride (MgB2) recently made an impact as a 
Jpjromising new superconductor with a surprisingly high crit- 
Qal temperature Tc for a simple binary compound. Besides 
temperature T, the other principal parameters that define 
the operational space of a superconductor are magnetic field 
Wif .and current density j. Superconductivity perishes above 
i^c2(T, j) and jd{T,H). For most practical applications it is 
Jw)t sufficient for the system to be merely in a superconduct- 
^f^yg state (presence of a finite order parameter amplitude) but 
fr^ r the system to be in a dissipationless state (constancy of 
der-parameter phase difference across sample length) . Thus 
C^om a practical standpoint, the conventional critical current 
^3fensity jc is limited by the depinning of vortices and their 
+6eiisequent motion. In the high- Tc cuprates, the large Tc is 
[:companied by high values^ of jd (> 10® A/cm^). However 
(le layered structure and small coherence length ^ lead to very 
Sjeak vortex pinning. Thus between jc and jd, there can be a 
road dissipative regime?:. 
Q MgB2 has an intermediate Tc and low values of jd and Hc2- 
^ ~ 10^ A/cm^ at low temperatures and, in single crystals, 
j^c2 ~ 3 T parallel to the c-axi&*,^ the films studied here, the 
• Trttrinsic disorder makes two striking improvements. First of all 
^^c2 is enhanced six-fold. Secondly, pinning is enhanced to the 
libint of raising jc to the same order of magnitude as jd- These 
■ observations are consistent with the reduction of ^ with disor- 
der. Finally MgB2's normal-state resistivity p„ is lower than 
that of a typical cuprate superconductor (e.g., YiBa2Cu307) 
by two orders of magnitude. This increases the vortex vis- 
cosity by the same factor, so that even when the vortices are 
depinned, they flow at very low velocities causing insignificant 
dissipation until jc becomes comparable to jd- AH these fac- 
tors lead to an extremely steep non-linear current- voltage {IV) 
curve with a complete absence of an Ohmic regime character- 
istic of free flux flow^. 



II. EXPERIMENTAL DETAILS 



The samples are 400 nm thick films of MgB2 fabricated using 



An amorphous boron film was deposited on a (1102) AI2O3 
substrate at room temperature by pulsed-laser ablation. The 
boron film was put into a Nb tube with high-purity Mg metal 
(99.9%) and the Nb tube was then sealed using an arc fur- 
nace in an argon atmosphere. Finally, the tube was heated 
to 900° C for 30 min. in an evacuated quartz ampoule sealed 
under high vacuum. X-ray diffraction indicates a highly c- 
axis-oriented crystal structure normal to the substrate surface 
with no impurity phases. The films were photolithographically 
patterned down to narrow bridges. In this paper we show data 
on three bridges, labelled S, M, and L (for small, medium, and 
large) with lateral dimensions 2.8 x 33, 3.0 x 61, and 9.7 x 
172 /im^ respectively. The lateral dimensions are uncertain by 
±0.7/im and the thickness by ±50 nm. Fig. ^a) shows the 
sample geometry. The horizontal sections of the current leads 
add a ~ 15 % series resistance to the resistance of the actual 
bridge, which enters in measurements made at very small j. 
At high J in the mixed state, the extra contribution is frozen 
out because the current is spread out in these wider areas, di- 
luting its density, so that one observes mainly the resistance of 
the bridge (A more intricate contact geometry was prohibited 
because of the extreme difficulty we faced in etching MgB2.). 

The non-linear electrical transport measurements were made 
using a pulsed signal source with pulse durations ranging from 
0.1 to 4 fis and a duty cycle of about 1 ppm. At the other 
extreme, a conventional continuous DC method at a very low 
current (/ = 1.4/iA) was employed for the resistive traces used 
to determine iJc2- Fig-^b) shows pulse waveforms under the 
especially severe conditions of j = 9.7 MA/cm^, E = 83 V/cm, 
and p ^ jE ^ 803 MW/cm^. The resistance rises to 90% of 
its final value in about 50 ns from the 10% onset of /. From 
a knowledge of the thermal conductivities and specific heat 
capacities of the film and substrate materials, and their mu- 
tual thermal boundary resistance, one can calculate the total 
thermal resistance Rth for any pulse duration^. Also if R{T) 
has enough variation, the film's own resistance can be used 
as a thermometer to measure Rth- For films of YiBa2Cu307 
(YBCO) on LaAlOs, which were used for most of our previ- 
ous work, we found Rth ~ 1-10 nK.cm'^/W at microsecond 
timescale32i2ii£. In the case of our MgB2 films, we expect 
Rth to be smaller because of sapphire's very high conductiv- 
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FIG. 1: (a) Sample geometry used for resistance measurement. At 
low values of j the wide lead areas add a constant resistance of about 
15% of the total value. At high j this contribution is frozen out. (b) 
Pulse waveforms under worst-case conditions (j — 9.7 MA/crr?', 
£; = 83 V/cm, and p = jE = 803 MW/cm^ on the plateaus). The 
resistance rises to (90% of) its final value in about 50 ns from the 
(10%) onset of I. 



first principles are not all known for this film-substrate combi- 
nation and MgB2 has a very flat R{T) below 50 K, so one can't 
measure Rth as was done for YBCO. We can, however, obtain 
an upper bound on Rth in the following way: Fig. EJa) shows 
IV curves for sample L in zero field (This is the largest sample 
with the lowest surface-to-volume ratio, so that it represents 
the worst case thermal scenario.). The curves were measured 
with the sample in different thermal environments. Above 
some threshold current Id ~ 650 mA, the system abruptly 
switches into the normal state. The value of Id is not sensitive 
to the thermal environment contacting the exposed surface of 
the film, confirming that the highly conductive sapphire, to- 
gether with the greatly reduced heat input during the short 
pulse, prevents the film's temperature from rising significantly 
(It has been shown by StoU et ahW that if there is sample 
heating, the thermal environment makes a significant differ- 
ence because it will provide an additional path for the heat to 
flow through.). Fig.|21[b) shows the top portion of one of the IV 
curves. The resistance jumps directly to the full normal-state 
value (The arrow indicates the first data point with non-zero 
resistance; the dashed line corresponds io V ~ Rnl )■ It is 
argued elsewherei^ that this jump to the normal state occurs 
due to pairbreaking by the currenti^. At the point the system 
is just driven normal, the power density reaches p = jE = 1.01 
GW/cm^ (arrow in Fig. |2Ib)). This sets a gross uppcrbound 
of Rth ~ 7 nK.cm'^/W. Note that the main bottle neck of heat 
conduction is the film-substrate boundary resistance which is 
not strongly temperature dependent®. In the present work, 
typical p values are two orders of magnitude lower than the 
critical 1 GW/cm'^ and so we expect the temperature rise to 
be a small fraction ('^l %) of Tc. 

The magnetic field is applied normal to the film (parallel 
to the c axis) and the self field of the current (< 50 G) is 
much lower than the applied fields used in this work. Further 
details of the measurement techniques have been published in 
a previous review articlw^ and other recent papergSiiS. 



III. RESULTS AND ANALYSIS 

Fig. O shows the resistive transitions at a low continuous 
current of / = 1.4^A in different fixed magnetic fields. Panel 
(a) shows the full curves for sample M, whereas for sample S we 
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FIG. 2: Zero-field IV curves for sample L show an abrupt switch 
to a state of full normal-state resistance, R„, by currents of pair- 
breaking magnitudes, (a) The switch to a dissipative state is not 
significantly influenced by the thermal environment, (b ) Top portion 
of the curve in superfluid He^ . The resistive portion of the measured 
curve extrapolates to the V = Rnl dashed line. 



full and central views look very similar for all three samples. 
In the central region, the shifts are approximately parallel, 
reducing the ambiguity of Hc2(T) defined by the transition- 
midpoint criterion. Fig. 0] shows those IIc2{T) values. IIc2{T) 
has an unusual linear dependence that departs noticeably from 
the dashed line of the WHH (Werthamer-Helfand-Hohenberg) 
functiorJS,. This observation is consistent with that of a dirty 
two-gap superconductor—. Panel (b) shows Hc2 {T) for sample 
S defined by other criteria (see Fig.lSJb) above). Note that the 
linear dependence is preserved and the slopes (and hence the 
inferred zero-T values of iJc2 ) do not change significantly with 
the choice of criteria. We infer a value of 17± 3 T for iJc2(0), 
which is six times higher than the value of 3 T found in single 
crystals'*. This can be understood in view of the much higher 
resistivity and hence shorter coherence length ^: For a dirty 
superconductor ^ « VS.ol^ where is the clean-limit value of 
^ and I is the impurity scattering mean-free-path. Since p oc I 
and Hc2 oc 1/^^, we expect i/c2 oc p. The normal resistivity 
of our film samples (14 /Ltf2-cm) is about seven times higher 
than the value (2 ^Jl-cm) of the single crystals of Refi^. Note 
that this Hc2 enhancement occurs naturally in the as-prepared 
films without any special effort to introduce pinning sites. The 
microscopic cause of the enhanced scattering leading to the 
higher i7c2 in these samples is presently not clear. 

We now turn to the in-field IV curves to investigate the na- 
ture of flux motion. In a system with weak flux pinning, the 
resistance goes through alternate regimes of Ohmic {V oc /) 
and non-Ohmic behavior. At very low driving forces (low j) 
there can be observable resistance due to thermally activated 
flux flow (TAFF) or flux creep. Then one encounters a non- 
linear response as current-driven depinning sets in; in effect 
the number of mobile vortices is rising. This is incipient flux 
flow. At sufficiently larger j, the vortex motion is effectively 
free from the influence of pinning and the response becomes 
Ohmic again. We previously introduced the term free flux 
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FIG. 3: Resistive transitions of two MgB2 bridges at indicated flux 
densities. 1=1.4 M^- (O') Sample M; midpoint Tc's lie on the dashed 
line through R = Rn/'2 = 1.95Q. (b) Similar set for sample S. 
Dashed lines run through R = Rn/2 = 1.22Q., and also through the 
lower and higher resistive criteria R = 0.840 and R = 1.76n. 



resistivity should correspond to the canonical Pf ^ PnB/Hc2 
Bardeen-Stephen expressioni^. As one goes beyond FFF, non- 
linearity can set in because of heating of the electron gas^ or 
changes in the electron distribution function^"*. Finally at yet 
higher currents, pair-breaking destroys superconductivity and 
drives the system normal. Here the resistance again ceases to 
change with current, being characteristic of the normal state, 
and so the response becomes Ohmic one more time. These 
stages of dissipation have been described in our previous re- 
view article. In YiBa2Cu307, the depinning critical current 
is sufficiently weak compared to the pair-breaking value that 
all of the regimes can be observed. 

In MgB2 the situation is very different. Fig. E^a) shows the 
R(I) curves of the present MgB2 samples. After the onset of 
dissipation, the resistance quickly rises to the full normal-state 
value. It should be noted that the plateaus do not correspond 
to FFF but to the normal state. Accordingly the resistance 
value changes very little with the applied B, especially for the 
curves at lower fieldsii. The overall shapes of the curves are 
almost independent of field. When the curves of the previ- 
ous panels are shifted vertically and horizontally by constant 
amounts they can be made to overlap as shown in Fig.|SJb). 

Fig. IHl shows the raw and shifted curves for the other two 
samples. 

These steep IV characteristics in MgB2 are similar to those 
for niobium alloys but are markedly different from those in the 
cuprates. The latter can be understood in terms of MgB2's 
lower temperature scale and much greater pinning (due to 
its higher isotropy and ten times larger vortex cross section). 
Hence thermal activation and current driven depinning are de- 
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FIG. 4: Hc2 versus B for three samples. The solid straight lines 
are linear fits to the data. The dashed lines correspond to the WHH 
function. Panel (b) shows Hc2 for sample S defined at two other 
criteria. The linear behavior of Hc2{T) continues to hold for the 
R=0.84 ^ and R=1.76 SI criteria. 
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FIG. 5: Resistance versus current vurves for sample M. Flux den- 
sities are indicated from left to right. The sample was immersed in 
superfluid helium and T = 1.5 K for all data, (a) Raw data, (b) 
Same data linearly shifted so as to collapse onto a single common 
curve. 



response to become steeply non-linear, rising from very lit- 
tle dissipation to the full normal state within a rather narrow 
range of currents. The collapse in Figs. Eland |S1 seems to sug- 
gest that the rise in resistance stems mainly from the intrinsic 
mechanisms of current- and field-induced pair-breaking, rather 
than from current driven depinning, even when B approaches 
Hc.2- This makes MgB2's mixed-state response unique among 
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FIG. 6: Resistance versus currenv curves for samples S and L. Flux 
densities are indicated from left to right. The sample was immersed 
in superfiuid helium and T — 1.5 K for all data. Panels in the left 
column show raw data and panels in the right column show linearly 
shifted data. 



IV. CONCLUSIONS 



We have investigated the low-temperature {T <C Tc) in- 
field transport behavior of MgB2 and present the first mea- 
surement of the full dissipation curves (i.e., < j jd and 
< Rt=o ^ Rn) for this system. MgB2 films made by the 
two-step laser-ablation process have an intrinsic pinning of a 
critical magnitude, such that the principle dissipation and re- 
sulting IV curves arise mainly from intrinsic mechanisms such 
as pair-breaking. The onset of dissipation is within an order 
of magnitude of the pair-breaking current, even at flux den- 
sities of a few teslas — the resistance rising quickly to the full 
normal-state value as the current is increased beyond jc- The 
same disorder that enhances pinning, also enhances i?c2 and 
qualitatively changes its temperature dependence because of 
the two-band nature of the superconductivity in this material. 
As explained by the theory of Gurevich, Hc2{0) extrapolates 
to a higher value than would be expected for WHH behav- 
ior, and the slope, dHc2/dt, is relatively constant. Our mea- 
surements of Hc2 are consistent with this predicted unusual 
dependence. Both the enhancement in Hc2 and the critically 
pinned flux make such MgB2 films more promising for applica- 
tions, besides providing a cleaner and more intrinsic view of a 
superconductor's current-voltage response in the mixed state. 
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The value of the current density at which this jump occurs, ~2 x 

10^ A/cm^, is roughly comparable to the theoretical estimate of 

jd(0) = cHc(0)/[3V67rA(0)] ~ 6 X 10^ A/cm^. 

Sometimes large departures can occur for exceptional situations 

such as superclean systems and narrow vortex cores where the 

internal energy-level spacing exceeds their widths. 

The slight shift in plateau resistance at the highest fields can 

be understood in terms of spreading of resistance outside the 

bridge area and into the current-lead areas as explained in the 

experimental section. Fields approaching Hc2 start driving the 

whole film normal at relatively low currents so that the resistance 

of the wider current-lead areas is not frozen out. This causes the 

normal-state plateau to rise slightly at the highest fields. 



